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Abstract 

We derive explicit formulas for calculating e A , cosh A, sinh^, cos A 
and sin A for a given 2x2 matrix A. We also derive explicit formulas 
for e A for a given 3x3 matrix A. These formulas are expressed 
exclusively in terms of the characteristic roots of A and involve neither 
the eigenvectors of A, nor the transition matrix associated with a 
particular canonical basis. We believe that our method has advantages 
(especially if applied by non-mathematicians or students) over the 
more conventional methods based on the choice of canonical bases. We 
support this point with several examples for solving first order linear 
systems of ordinary differential equations with constant coefficients. 
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1 Introduction 



The exponential e A of a square matrix A and the related one-parameter fam- 
ily e tA are important concepts in mathematics. Here is one example (among 
many others): Let x'(t) = Ax(t) be a system of first order homogeneous or- 
dinary differential equations with constant coefficients with initial conditions 
x(0) = xq G M. n , where A is an n x n matrix with real entries. Then the so- 
lution to the system is given by the formula x{t) = e tA xq (Michael Artin [1], 
p. 140). 

The exponential e tA is defined by the Taylor expansion: e tA = Y2^=o ^TTp- 
For some particular matrices A the exponential can be easily calculated. For 
example, if N is a nilpotent matrix of order m, i.e. N m = O, then e tN = 

En-i ^ If p is a projectioil) Le . P 2 = Pt then *P = j + *r = 

I—P+e t P. If A = al for some scalar a, then e tA = Y^°° n 1 = e at L From 

Z — Jn=U 77,; 

the Taylor expansion it follows that if AB = BA for two square matrices, 
then e A+B = e A e B . Also, if A and C are similar matrices with a transition 
matrix T, i.e. A = TC T" 1 , then e A = Te c T^ 1 . 

These facts can be found in many textbooks (Jerry Farlow et al [2], p. 
350) and handbooks p. 132-133) on linear algebra. For a general ma- 
trix A, however, the rule for calculating e tA becomes somewhat more com- 
plicated (A.I. Malcev [3], p. 118 or Michael Artin P, p. 480-482): We 
have to find the Jordan form T _1 AT = diag(Ci, C2, ■ ■ ■ , Ck) of the matrix 
A, where Cj's are Jordan cells. Thus the formula e A = Te c T~ 1 becomes 
e A = Tdiag(e Cl , e° 2 , . . . , e Ck ) T~ l . On the other hand, each cell can be pre- 
sented in the form Ci = \I + N, where N mi = O, Aj is the corresponding 
eigenvalue of A and m ; = size(Cj). As a result, e Ci = e Xi Y^L® 1 ^T- The 
procedure is universal and it looks attractive but it has the following disad- 
vantages: (a) The algorithm for calculating the transition matrix T and its 
inverse T _1 is difficult and time consuming. It is especially complicated in 
the case of multiple characteristic roots, when the calculations of the general- 
ized eigenvectors require skills that are too advanced for students in a typical 
course on ordinary differential equations and linear algebra, (b) The origi- 
nal framework R n must be extended to C n , which again might be confusing 
for students. Our observation is that students and non-mathematicians (the 
latter might have rusty knowledge of linear algebra) who rarely calculate ex- 
ponentials of matrices often have considerable difficulty calculating e tA even 
for 2x2 and 3x3 matrices. 

In this article we derive explicit formulas for calculating e tA for a 2 x 2 
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and 3x3 matrix A with real entries. In the case of 2 x 2 matrices A we also 
derive explicit formulas for cosh A, sinhA,cosA and sin A. These formulas 
are expressed exclusively in terms of the characteristic roots of A and in- 
volve neither the eigenvectors of A, nor the transition matrix associated with 
a particular canonical basis. We believe that our formulas are suitable for 
handbooks in the sense that one can apply them with limited background in 
linear algebra, in particular, without having slightest idea what a "transition 
matrix is". The formulas are derived with the help the characteristic poly- 
nomial Ja (A) = det(A — XI) of A and the Cayley-Hamilton theorem which 
says that /a(^4) = 0. But, again, the formulas can be applied without knowl- 
edge of the Cayley-Hamilton theorem. In Corollary 12.11 and Corollary 14. II we 
identify some matrices closely related to A which allow easy calculation of 
the invariant subspaces of A and the canonical bases (if needed). We rec- 
ommend our method for teaching exponentials in a course on linear algebra 
and ordinary differential equations. To support this point we present several 
examples for solving linear systems of ordinary differential equations. 

For other methods for calculating the exponential of a matrix not men- 
tioned in our paper we refer to the survey article (Cleve Moler and Charles 
Van Loan [1]), where the reader will find more references to the subject. 

In what follows / and O denote the identity and zero square matrices, 
respectively. We denote by [x\y \z \ ...] the matrix with column- vectors 
x, y , z, • ■ ■ G M n . Also, E\ = ^(A) denotes the eigenspace of A corresponding 
to ju 6 C. 

2 Exponential of a 2 x 2 Matrix 

We derive formulas for calculating the exponential of a given 2 x 2-matrix A 
based exclusively on the characteristic roots of A. Our formulas are easy to 
memorize and simple to use in the sense that: (a) our framework is always IR 2 
(never C 2 ); (b) our formulas involve linear operations between matrices only; 
(c) our formulas can be handled with limited background in linear algebra - 
they involve neither the eigenvectors of A, nor a particular canonical basis 
(although we need the characteristic roots of A). We illustrate this point 
with several examples of linear systems of ordinary differential equations. 
Although we do not need a change of the basis to calculate the exponential 
e tA , we also derive explicit formulas for a canonical basis and a transition 
matrix T for which C = T~ X AT is exceptionally simple for the purpose of 
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the exponentiation (Corollary 12. ip . 



Theorem 2.1 Let A be a 2 x 2-matrix with real entries and let Ai and A 2 
be the characteristic roots of A. 

Case 1: If Ai = A 2 (real), then A = X I + N, where A = Ai = A 2 and 
the matrix N = A — X Q I satisfies the equation N 2 = O. Consequently, for 
every real (or complex) t we have 

(1) e tA = e Xat (I + tN). 

Case 2: If Ai )2 = ct ± iuo for some a, u G R, u 7^ 0, t/ien A = ai + uoJ , 
where the matrix J = -(A — ai) satisfies the equation J 2 = —I. Conse- 
quently, for every real (or complex) t we have 

(2) e tA = e at [(cos ut)I + (smut) J}. 

Case 3: If Xi ^ X 2 (both real), then A = aI + (3J, where a = f3 = 

Al ~ A2 and the matrix J = \{A — ai) satisfies the equation J 2 = I . Conse- 
quently, for every real (or complex) t we have 

(3) e tA = e at { (cosh/%)/ + (sinh (3t)J } = X - [e Xl \l + J) + e X2t (I - J)] . 

Proof: Case 1: We first show that N 2 = 0. Indeed, N 2 = (A - X I) 2 = 
A 2 — 2XqA + Aq/ = O, by the Cayley-Hamilton theorem, since 2Ao = tr(A) 
and Aq = det(A). Thus we have 

e tA = e t(x i+N) = e x tl . e tN = e A t 7 g + tN + + ...) = e Xot (I + tN), 
as required. 

Case 2: We first show that J 2 = —I. Indeed, we have A 2 — 2al + (a 2 + 
uj 2 )I = O, by the Cayley-Hamilton theorem, since 2a = X\ + A 2 = ti(A) and 
a 2 + u 2 = AiA 2 = det(A). Thus J 2 = (~{A - ai)) 2 = ^(A 2 - 2al + a 2 I) = 
4r (A 2 - 2 7 J + (a 2 + uo 2 )I - to 2 1) = MO- uo 2 I) = —I. Next, we calculate 



J 00 / . + \2n j2n 00 ( .4-\2n+l j2n+l 

(2n)\ ^ (2n+l)\ 

n=0 v ' n=0 v ' 



t A 



(2n)! ^ 

n=0 v 7 n=0 

e ™^ + e ™T4 - + (^), } 



(2n)! (2n + l)! 



n=0 v ' n=0 
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as required. 

Case 3: We shall, first, show that J 2 = /. Indeed, we have A 2 — 2al+(a 2 — 
f3 2 )I = O, by the Cayley-Hamilton theorem, since 2a = X\ + A2 = tr(A) and 
a 2 - (3 2 = AiA 2 = det(A). Thus J 2 = - al)) 2 = jp(A 2 - 2al + a 2 I) = 
jz(A 2 - 2al + (a 2 - (3 2 )I + (3 2 I) = 4(0 + (3 2 I) = I. Next, we calculate 



f 00 ( o + \2n j2n 00 ( o f \2n+l T 2n+1 
t(aI+(3J) = p atl p f3tJ = p at ) l^j J , l/ 3 ^ J 



//I 



(2n)! ^ (2n+T 

n=0 v 7 n=0 v ' 



E W + E = e " {(cosh/3t)J + (sinh/%) J} 



(2n)\ ^(2n + l 

n=0 v 7 n=0 v 



as required. ▲ 

The formulas ([I])-© show that we do not need the eigenvectors of A in 
order to calculate e tA . However, the eigenvectors of A as well as canonical 
forms and canonical bases (if needed for classification, graphing, etc.) can 
be easily extracted from the matrices N and J. The next corollary follows 
easily from the above theorem and we leave the proof to the reader. 

Corollary 2.1 (Eigenvectors, Canonical Forms and Bases) Under the 
notation of the previous theorem we have the following: 

Case 1: If N ^ O, then E\ = \ (A) = Im(iV). In particular, either of the 

Ao 1 \ rp-l 



non-zero columns ofN is an eigenvector of A. We have A = T ( ^ ^ j T 

where T = [Nx\x] and fel 2 such that Nx ^6. If N = 0, then A = X I. 
Case 2: The matrix A does not have eigenvectors in R 2 . We have 

A = T ( a ] T _1 , where T = \ Jx I x] and x G R 2 , x ^ 0. Notice that 
\—lo a J iii 

the matrix ( a ^ is conformal (i.e. the corresponding transformation 

preserves the angles in M. 2 ). Alternatively, in C 2 the matrix A is diagonaliz- 
able with eigenvalues a ±iu and eigenvectors (I =F i J)x, f / 0, respectively. 
Case 3: We have E X= \ 12 (A) = Im(J ± J), respectively. We have A = 

T ® j T~ x , where T = [(I + J)x \ (I — J)y) with x, y e R 2 such that 

(I + J)x^6 and (I - J)y ^ 0. 

Here are several examples tested in class by the second author in a course 
on linear algebra and ordinary differential equations. We strongly recommend 
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these formulas for teaching exponentials. 

Example 2.1 We shall find the solution of the initial value problem without 
involving eigenvectors: 

x' = 3x + 2y, 
y ' = -8x - by, 

with initial conditions x(0) = 1, y(0) = —1. We have A = ( ^ ^ 
and xq = ( J . The characteristic polynomial is A 2 + 2A + 1 wrf/i roots 



x -1 / 

Ai = A 2 = —1 (Case 1). We calculate 



N = A + I 



3 2 \ A 0\ _ / 4 2 
-8 -5 J + VO 1 7 V —8 -4 



We leave to the reader to verify that N 2 = O. For the exponential of A we 
apply formula (TTJj: 

77ms x = e~*(l + 2t), y = -e"*(l + 4t). 

Remark 2.1 //a canonical form for A and a transition matrix are still 
needed, we can easily calculate them with the help of Corollary \2.1[ For 

example, T~ X AT = f ^ ^l)' w ^ ere ^ = — ( ^ f or % — 

(0, 1). Notice that either of the columns of N is an eigenvector of A. 

Example 2.2 We shall find the solution to the initial value problem without 
involving eigenvectors: 

x = y, 

V = -5ac - 2y, 

with initial conditions x(0) = 2, y(0) = 1. We have A = ( ^ and 
Xq = ( ? J ■ 27ie characteristic polynomial is A 2 + 2A + 5 wift roofo Ai^ = 



-1 ± 2%, thus a = — 1 and uo = 2 (Case 2). We calculate 



J 



■(A-al) 



1 

-5 -2 



1 
1 



1 / 1 



We leave to the reader to verify that J 2 = —I . For the solution to the system 
we apply formula (0|): 



-> tA-> 

x = e xq 



cos 



cos 2t + | sin 2t, 
-|sin2t, 



\ sin 2t 
cos 2t — \ sin 2t 



_ t ( 2 cos 2t + | sin 2t 
cos 2t - y sin 2t 



Thus x = e _t (2 cos 2t + | sin 2t ) , y = e - * (cos 2t - ^ sin 2t) . 

Remark 2.2 //a canonical form for A and the corresponding transition 
matrix are still needed, we can easily calculate them with the help of Corol- 

1 2 N 
-5 0. 


-1 - 2il ' 



laru \2.1\ For example, T 1 AT 



2 -y 1 , where T = [Jx \ x 



forx = (2,0). In the framework of 'C 2 we have T l AT 



■l + 2i 




where T = [{I — iJ)x \ (I + % J)x 



2 + i 2-i 
— 5i 5i 



forx = (2,0). 



Example 2.3 We shall find the solution to the initial value problem without 
involving eigenvectors: 

x' = 5x — y, 
y' = 3x + y, 



with initial conditions x(0) = 1, y(0) = 2. We have A 



5 -1 
3 1 



and Xq 



. The characteristic polynomial is A 2 — 6A + 8 with roots Ai = 4, A2 = 2 



(Case 3). We calculate a = 3, j3 = 1 and 



J=~(A-aI) 



5 -1 
3 1 



3 
3 



2 -1 

3 -2 
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We leave to the reader to verify that J 2 = I . For the solution to the system 
we apply formula 

( cosht + 2 sinht, — sinht \ (\ 



x(t) e tA [ 2 J e 3t y 3 ginh ^ cogh t _ 2 ginh t j y 2 



e 3t 



cosht 
2 cosh t — sinh t 



Thusx = e 3i cosht = ±(e 4 ' + e 2t ) and y = e 3t (2cosht-sinht) = |(e 4 ' + 3e 2 '). 

Remark 2.3 If a canonical form for A and the transition matrix are still 
needed, we can easily calculate them with the help of Corollary \2.1[ For 

example, T~ l AT = 2 J ' w ^ ere T — [(I + J)x \ (I — J)x] = ( j /or 

x = (0, 1). Notice that the columns ofT are eigenvectors of A. 

3 Some Applications 

If A is a square matrix, we define the matrices cosh A, sinh iA, e lA , cos A and 
sin A by the corresponding Taylor series. For example, cosh A = J2^=o finYp 
and the rest are defined similarly. In this section we derive explicit formulas 
for these elementary functions which follow easily from Theorem 12.11 

Corollary 3.1 (Hyperbolic Functions) Let A be a 2 x 1-matrix with real 
entries and let Ai and A2 be the characteristic roots of A. 
Case 1: If X\ = X2 (real), then 

cosh A = (cosh X )I + (sinh Xq)N, sinh A = (sinh X )I + (cosh Xo)N, 

where Xq = Ai = A2 and N = A — XqI. 

Case 2: If X\^ = a ±iu for some a, uo G K, c<j 7^ 0, i/ien 

cosh A = (cosh a) (cos a;)/ + (sinh a) (sin a;) J, 
sinh /I = (sinh a) (cos cj)/ + (cosh a) (sin u>) J, 

where J = -(A — ai). 

Case 3: If X\ 7^ A2 (both real), then 

cosh A = (cosh a) (cosh/?)/ + (sinh a) (sinh (3) J, 
sinh A = (sinh a) (cosh (3)1 + (cosh a) (sinh (3) J, 

y,/iere a = ^±^, /3 = and J = ±(4 - a/). 
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Proof: The results follow directly from Theorem 12 . 1 1 and the formulas cosh x = 
\{e x + e~ x ) and sinhx = \{e x — e~ x ). A 

Corollary 3.2 (Complex Exponent) Let A be a 2 x 2-matrix with real 
entries and let Ai and X 2 be the characteristic roots of A. 
Case 1: If X\ = X 2 (real), then 

e ±iA = e ±iAo( /±?7V ) ; 

where Xq = Ai = A2 and N = A — XqI. 

Case 2: If Xi >2 = a ±iu for some a, uo G K, uo 7^ 0, then 

e ±iA = e ±i«[( cosha; )/ ± i( s inhcj) J], 

where J = -(A — oil). 

Case 3: If X\ 7^ A2 f&oi/i rea/j, t/ien 

where a = ^±^> /3 = anrf J = U A ~ aI )- 

Proof: These formulas follow directly from Theorem 12.11 for t = i. A 

Corollary 3.3 (Trigonometric Functions) Let A be a 2 x 2-matrix with 
real entries and let X\ and X 2 be the characteristic roots of A. 
Case 1: If Ai = A2 (real), then 

cos A = (cos Xo)I — (sin Ao) N, sin A = (sin X )I + (cos Ao)iV, 

where Xq = Ai = A2 and N = A — XqI. 

Case 2: If X\ i2 = a ±iu for some a, uo G K, c<j 7^ 0, i/ien 

cos A = (cosa)(coshc<j)/ — (sin a) (sinh o>) J, 
sin A = (sin a) (cosh u>) I + (cosa)(sinhu;) J, 

where J = -(A — ai). 

Case 3: If X\ 7^ A2 (both real), then 

cos A = (cos a) (cos/?) / — (sin a) (sin/?) J, 
sin A = (sin a) (cos/?)/ + (cosa)(sin/3) J. 

where a = ^±**, /3 = and J = ±(4 - a/). 

Proof: As before, the result follows directly from Corollary 13.21 and the for- 
mulas cos a; = |(e lx + e~ tx ) and sin a; = ^-{e lx — e~ lx ). A 
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4 Exponential of a 3 x 3 Matrix 



In this section we extend the result about e A , presented in Section [2, for 
3x3 matrices A. We should notice that: (a) the hardest (and the most 
time consuming) parts of our formulas for e A are the calculating the squares 
of some matrices; we need neither to solve linear systems of equations, nor 
to invert matrices: (b) our framework is always R 3 (never C 3 ). For that 
reason we believe that our method has advantages over the most conventional 
methods based on the choice of a particular canonical basis with transition 
matrix T and calculating its inverse T _1 . To support this point we present 
several examples for solving linear systems of ordinary differential equations 
in three unknowns. 

In what follows we shall use the fact that if A is a matrix and B = A — ai 
for some real a, then /s(A) = /a(A + a). 

Theorem 4.1 Let A be a 3 x 3 matrix with real entries and let Ai,A2 and 
A3 be the characteristic roots of A. Then: 

Case 1: If Ai = A2 = A3 = Ao, then N = A — X I is a nilpotent matrix 
such that N 3 = 0. Consequently, for every t el we have 

(4) e tA = e Xot (l + tN + ^N 2 

Case 2: If Ai = A 2 ^ A 3; then A = X I + (A 3 - A )P + iV, where 
Ao — Ai = A2, 

' (A - A /) 2 and N = A — A / - (A3 - X )P. 



(A 3 - A ) 



2 



Also, we have PN = NP = 0, P 2 = P and N 2 = 0. Consequently, for 
every t e 1R we have 

(5) e tA = e Xot (I - P + tN) +e Xst P. 

Case 3: Let Ai,2 = ct ± iui for some a,w 6 R, w 7^ 0. Then A = 
ai + (A 3 — a)P + ijj J, where 

(A - ai) 2 + u 2 I I 

P = ^~\ ^2~i — 2~> J = - [A-al - {X 3 -a)P}. 

(A 3 — a) 2 + u> 2 uj 
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Also, we have PJ = JP = 0, P 2 = P, J 2n = (-l) n (/ - P), n = 
1,2,..., J 2n+1 = (-l) n J, n = 0,1,2,.... Consequently, for every t E R 
we have 

(6) e tA = e at [ (cos ut) (I - P) + (sin ut)J] + e Xat P, 

Case 4-' If^i, A2, A3 are three distinct reals, then A = aI+(Xs—a)P+(3J, 
where a = (Ai + A2V2, (3 — (Ai — A2V2, and 

P = Tx ^ = ^ - «/ - As - a P . 

Abo, PJ = JP = 0, P 2 = P, J 2n = I - P, n = 1, 2 . . . , J 2n+1 = J, n = 
0,1,2.... Consequently, for every t EM. we have 



(7) e iA =e at [(cosh /ft) (J - P) + (smh(3t)J\ + e As *P = 

^ I e Ait(/ + J - P) 4- i e A2t (7 - J - P) + e A3 *P. 

Proof: Case 1: We have f A (X) = (A - A ) 3 and iV 3 = (A - A /) 3 = f A {A) = 
0, by the Cayley-Hamilton theorem. Next, e tA = e * (A ° /+Ar) = e tXo1 e tN = 
e A( '*/(J + tJV + (t 2 /2)iV 2 + + . . . ) = e Xot (I + tN+ (t 2 /2)N 2 ), as required. 

Case 2: Denote A 3 — A = b, A — X I = B and observe that P = ^P 2 and 
N = B-bP = -\B{B-bI). Notice that f A {\) = (A-A ) 2 (A-A 3 ) implying, 
by translation, / B (A) = A 2 (A - 6). It follows P 2 (P - bl) = f B (B) = 
(by the Cayley-Hamilton theorem), i.e. P 3 = 6P 2 and P 4 = b 2 B 2 . After 
these preliminary evaluations we calculate P 2 = P 4 /6 4 = b 2 B 2 /b 4 = P, as 
required. Next, we have NP = PN = -^P 2 P(P - bl) = -^Bf B (B) = 0, 
and iV 2 = -^P 2 (P - bl) 2 = ±f B (B)(B - bl) = 0, as required. Finally, 

e tA = e t[X I+bP+N] = e X tI e btP e tN = 
t n b r, 



e x tj I j + P J^—j (/ + tiV ) = e ^ot[ 7 _ P + e 6*p](/ + t N) = 
e \ t(j _ p + t N) + e Xst P, 



as required. 

Case 3: We denote A — ai = B and A3 — a = b and observe that P = 
^(B 2 + u 2 I) and J = ±(B — bP). We have f A ( A) = [(A-a) 2 + cu 2 ](A- A 3 ) 
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which implies /s(A) = (A 2 + cu 2 )(A — b). By the Cayley-Hamilton theorem, it 
follows (P 2 + to 2 I)(B - bl) = implying also P 3 + co 2 B = b{B 2 + cu 2 /) and 
B 4 +u 2 B 2 = b 2 (B 2 +u 2 I). Next we calculate PB = BP = w ^ s (B 3 +u 2 B) = 
(P 2 + to 2 1) = bP. It follows B 2 P = b 2 P which helps us to show that 
P is a projection, i.e. P 2 = p^(P 2 + u 2 I)P = ^(B 2 P + u 2 P) = 
j^i^P + ^P) =P and also PJ = JP = ±(B-bP)P= ±(bP-bP) = 0. 
Finally, we have J 2 = ^(B-bP) 2 = ^(B 2 -2bBP + b 2 P) = ^(B 2 -b 2 P) = 
■^i[(b 2 +u 2 )P — u 2 I — b 2 P] = — (/ — P), as required. The rest of the formulas 
for the powers of J follow immediately. Next, we calculate e^ tJ : 



n=0 n=l v y n=0 v y 



V ; ^Z^ ( 2n ) v ~Z_> 2n + l 

n=l v ; n=0 v ; 

= P + (cos cut) (7 - P) + (sin cut) J. 
Finally, we calculate 

e tA = ^[al+bP+^J] = e atl e MP e vtJ = e at^j _ p) + & Up^ x 

x [P + (cos wt) (/ - P) + (smut) J]. 

The last leads to formula ([6]) after standard manipulations. 

Case 4: We denote, as before, A — al = B and A3 — a = b and observe 
that P = -^(B 2 - (3 2 I) and J = ~(B - bP). Next, we have f A (X) = 
(A — Ai)(A — A 2 )(A — A 3 ) which implies, by translation, /s(A) = (A 2 — /9 2 )(A — 6). 
It follows (P 2 — (3 2 I)(B — bl) = 0, by the Cayley-Hamilton theorem. Thus 
B 3 - (3 2 B = b(B 2 - (3 2 I) and B A - (3 2 B 2 = b 2 (B 2 - f3 2 I). Now the relations 
between P and J follow exactly as in Case 3. Next we calculate: 



e tA = e t[aI + bP + /3J] = e atl e btP ^tJ = e °*[J _ p + gWp] x 

[e at (/-P) + e X:it P]x 



X 



x [P + (cosh /ft) (J - P) + (sinh fit) J] = e at [(cosh /ft) (J - P) + (sinh fit) J] + 
e A3 *P, 

as required. ▲ 
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The formulas (JI])-(J7D show that we do not need the eigenvectors of A in 
order to calculate e tA . However, the eigenvectors of A as well as canonical 
forms and canonical bases (if needed for classification, graphing, etc.) can be 
easily extracted from the matrices N, P and J. The next corollary follows 
easily from the above theorem and we leave the proof to the reader. 

Corollary 4.1 (Eigenvectors, Canonical Forms and Bases) Under the 
notation of the previous theorem we have the following: 

Case 1: If N 2 ^ O, then E x=Xo {A) = lm(N 2 ). In particular, ei- 
ther of the non-zero columns of N 2 is an eigenvector of A. We have A = 

Ao i o\ 

T O A 1 T" 1 , where T = [N 2 x\Nx\x] and x G R 3 , N 2 x + 0. If 
\0 V 

N = O, then A = X I . If N 2 = O and N ^ O, then E x=Xo (A) ^ Im(iV). 
To complete a Jordan basis we can select an additional eigenvector by easy 
inspection of the columns of N. 

Case 2: IfN^O, then E X =\ (A) = Im(iV) and E X=X3 (A) = Im(P). In 
particular, either of the non-zero columns of N and P is an eigenvector of 

Ao i o\ 

A. We have A = T A 1 T" 1 , where 
\0 A 3 / 

T=[(A- X I)(A - X 3 I)x\(A - X 3 I)x\ Py], 

and x G R 3 , {A - X I){A - X 3 I)x + and Py + 0. If N = O, then 
A = diag(A , A , A 3 ). 

Case 3: The only eigenspace in R 3 is E X=X3 = Im(P). Here we have 
(a uo \ 

A = t\-uj a T~ x , where T = \J 2 x\Jx\ Py] and x,y eR 3 ,x ^ 0, 
^\ V 

Py 0. Alternatively, in C 3 the matrix A is diagonalizable with eigenvalues 
a ± iuj, A 3 and eigenvectors (J =F iJ 2 )x, f ^ 0, Py ^ 0, respectively. 

Case 4 : We have E x=Xl 2 (A) = Im(J±J 2 ), respectively, and E x=Xa (A) = 
Im(P). Thus the matrix A is diagonalizable in the basis consisting of the 
eigenvectors listed above. 

Here are several examples. 

2 -1 2 \ 

Example 4.1 Let x'{t) = Ax(t), where A = | 5 —3 3 , be a system 

-1 -2} 
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of three first order homogeneous ordinary differential equations with constant 
coefficients with initial conditions x(0) = (1, 0, 2). We apply the formula for 
the solution x(t) = e tA Xq. It remains to calculate e tA . We leave to the reader 



to check that Ai = A2 = A3 = — 1 (Case 1 in Theorem \4-l\ ). We calculate 

/ 3 -1 2 \ / 2 -1 1 \ 

N = A + 1 = 5 -2 3 , N 2 = ( 2 -1 1 \ , N 3 = O. 
V-l -lj V-2 1 -V 

The formula gives 

f ,2 n /l + 3t + t 2 , -t-t 2 /2, 2t + t 2 /2 

e tA = e \ t ( I + m+ _ N 2 \ = e -t 5t + t 2 ; 1 _ 2t — t 2 /2, 3* + t 2 /2 

V 2 / V -t-t 2 , t 2 /2, i- t -t 2 /2 / 

/1\ / 1 + 7t + 2t 2N 

T/ie solution of the system is x(t) = e tA I = e~* I lit + 2t 2 

\2/ \2 - 3t - 2t 2 y 

Remark 4.1 (Jordan Basis and Eigenvectors) We do not need a Jor- 
dan basis (and the eigenvectors) of A to calculate e tA . However, a Jordan 
basis for A, if needed, can be easily extracted from the above calculations with 
the help of Corollary \4.1\ For example, {N 2 e~z, Ne~3, e^} (where e~3 is the third 



vector in the standard basis in M. 3 ), forms a Jordan basis for A. Also, we 
have E\=-i = Im(A^ 2 ). In particular, either of the columns of N 2 is an eigen- 
vector of A for A = —1. We leave to the reader to check that A = TCT^ 1 , 

1-1 1 \ / 1 2 N 

where C = -1 1 J and T = [N 2 e" 3 \ Ne z \ e 3 ] = 1 3 

\0 -1/ V-i -1 1, 

The columns ofT forms a Jordan basis for A. 

(1 -3 

Example 4.2 Let x'(t) = Ax(t), where A — I 4 —7 8 | , be a system of 

\Q -7 7j 

three first order homogeneous ordinary differential equations with constant 
coefficients with initial conditions x(0) = (1,0, 1). We leave to the reader to 
check that Ai = A2 = — 1 = Aq and A3 = 3 (Case 2 in Theorem 4-l\ )- We 
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calculate 

, 2 



A + I 




A + I 
4 




N = A + I-4P= \ A -6 8-4 






-1 








1 






-2 




i-i 


-4 


2 






-2 


2y 




V-2 


1 





Notice that N 2 = O. The formula gwes 

i-2t t o\ /i -l r 

r / ! r ' | -4t l + 2t ] + (e 3 * — e _< ) [ 2 -2 2 
-2t t 1/ \2 -2 2, 

-2te-* + e 3 *, (1 + t)e _ * - e 3 *, -e~* + e 3t 

(-2 -4t)e"' + 2e 3 *, (3 + 2t)e- 4 - 2e 3 *, -2e~* + 2e 3 ' 
-4te"* + 2e 3 *, (2 + t)e _ * - 2e 3 ', -2e - * + 2e 3 ' 

/1\ /(-l-2t)e-* + 2e 3 * N 
For i/ie solution we have x(t) = e tA = (—4 — 4t)e _ * + 4e 3 * 

\1 / \(-2 - 4t)e~* + 4e 3i / 

Remark 4.2 (Jordan Basis and Eigenvectors) We just calculated e tA with- 
out the Jordan basis and eigenvectors of A. They, however (if needed), 
can be easily extracted from the above calculations with the help of Corol- 
lary \4-l\ For example, E\=-i = Im(iV) and E\=3 = Im(P). In particular, 
the first and second columns of N are eigenvectors of A for A = — 1 and 
either of the columns of P is an eigenvector of A and A = 3. We leave 

-1 1 N 

to the reader to check that A = TCT~ X , where C = I —10) and 





T=[(A-X I)(A-X 3 I)x\(A-X 3 I)x\Py] = ( 16 2 2 | forx = y = e x . 
The columns of T forms a Jordan basis for A. 

I 1 1 ~ l \ 

Example 4.3 Let x (t) = Ax(t), where A = I 3 \ , be a system of 

\1 l) 

three first order homogeneous ordinary differential equations with constant 
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coefficients with initial conditions x(6) = 
the standard basis in IR 3 . We have Ai >2 : 
(Case 3 in Theorem\4-l\). We calculate 



€2, where €2 is the second vector in 
- 1 ± i, i.e. 7 = to = 1, and A3 = 3 



P 



(A-^If + oo 2 ! 
(A 3 - 7 ) 2 + ^ 2 





J=-[A->yI-{X a ->Y)P] = Uo 
uj 

e Xst -e 7 * cos cut) P + e 7 * [(cos cot) 1 + (sin cut) J] 

5e*cost, 2e 3t — e*(2 cost — sint), — 5e*sint N 

0, 5e 3 *, 

5e*sint, e 3 * — e*(cost + 2 sint), 5e*cos£ 




For the solution we have x(t) = e tA 11 = k \ ^e 3 * 




2e 3 * — e*(2cost — sint) N 
5e 3 * 

e 3t _ e *( cos £ _|_ 2 sint) 

Remark 4.3 (Canonical Forms and Eigenvectors) As before, we just 
calculated e tA without the Jordan basis and eigenvectors of A. They, however 
(if needed), can be easily extracted from the above calculations with the help 
of Corollary \4-l In this case E\ = \ 3 = lm(P) and we have A = TCT^ 1 , 

(l 1 0\ /-I 2\ 

where C — I — 1 1 and T = [ J 2 x \ Jx\ Py] =10 5 for x = e[ 

\ 3/ \ 11/ 

and y = 5e"2. The columns ofT form a canonical basis for A. 
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